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Abstract 

We give a detailed and improved presentation of our recently pro- 
posed formalism for non-linear perturbations in cosmology, based on a 
covariant and fully non-perturbative approach. We work, in particu- 
lar, with a covector combining the gradients of the energy density and 
of the local number of e-folds to obtain a non-linear generalization of 
the familiar linear uniform density perturbation. We show that this 
covector obeys a remarkably simple conservation equation which is 
exact, fully non-linear and valid at all scales. We relate explicitly our 
approach to the coordinate-based formalisms for linear perturbations 
and for second-order perturbations. We also consider other quantities, 
which are conserved on sufficiently large scales for adiabatic pertur- 
bations, and discuss the issue of gauge invariance. 
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1 Introduction 



The relativistic theory of cosmological perturbations is an essential tool to 
analyze cosmological data such as the Cosmic Microwave Background (CMB) 
anisotropics, and thus to connect the scenarios of the early universe, such as 
inflation, to cosmological observations. Because the temperature anisotropics 
of the CMB are so small (5T/T ~ 10~ 5 ), considering only linear perturba- 
tions is an excellent approximation. This is why most of the efforts devoted to 
the theory of cosmological perturbations have dealt with linear perturbations 

013 EH HIS]. 

There are however some issues where one must take into account the non- 
linear aspects of cosmological perturbations. A good example is the study of 
inhomogeneities on scales much larger than the Hubble radius, where the Uni- 
verse could strongly deviate from a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) geometry. Another motivation, which has seen renewed interest re- 
cently, is to investigate the predictions of early universe models for primordial 
non-Gaussianity, with the hope to be able to detect this non-Gaussianity in 
the CMB data. This investigation requires the study of relativistic cosmo- 
logical perturbations beyond linear order |S10IB|- 

For small perturbations, as a first step beyond linear order, one can con- 
sider second-order relativistic perturbations. This is enough if one wants to 
compute the bispectrum of perturbations, an indicator of the presence of 
non-Gaussianities. In the context of general relativity, dealing with second- 
order perturbations is already an impressive task. The usual treatment 
introduces coordinates and the metric perturbations are then formally as- 
sumed to be the sum of a first-order quantity and of a second-order quantity 
P?l [Til frn [T^ [TTIl [T^l [T31 . One can treat similarly the fluid quantities and then 
write down Einstein's equations up to second order. As one can imagine, this 
approach is rather cumbersome. 

Other recent approaches of non-linear perturbations ^01 HZ1 (see also 
[19|) are based on the long wavelength approximation [2011^1122], which at 
lowest order is related to the so called separate universe picture that repre- 
sents our universe, on scales larger than the Hubble radius, as juxtaposed 
FLRW universes with slightly different scale factors j2Hl 121] ■ 

In the present work, we present the details of a different approach, which 
we proposed recently [2H], based on a purely geometrical description of the 
perturbations. Our approach is inspired by the so-called covariant formalism 
for cosmological perturbations introduced by Ellis and Bruni j2Hl- Although 
the linearized version of this formalism is often used and it has also been 
used at second order (see e.g. [211), our approach is fully non-perturbative 
and thus not limited to second-order perturbations. The equations and the 
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variables used in our approach encode the full non-linearity of cosmological 
perturbations on all scales. Although, at linear order, the covariant formalism 
is computationally equivalent to the much more used coordinate approach, 
the covariant approach turns out to be more efficient when one goes beyond 
linear order. 

In particular, we show that it is possible to define, in a geometric way, the 
generalizations of quantities widely used in the linear theory because they 
are conserved on large scales. These quantities are covectors which obey 
conservation equations with respect to the directional derivation (Lie deriva- 
tion) along the comoving worldlines. We give the full non-linear equations 
that govern these quantities. To emphasize the efficiency of our approach 
with respect to the more traditional coordinate-based approach, we show 
how gauge-invariant quantities, recently derived in the literature, can be ob- 
tained from our geometric quantities in a straightforward derivation. The 
particular case of scalar fields perturbations will be considered in a separate 
publication. 

This paper is organized as follows. In the next section, we give a brief 
overview of the covariant formalism for cosmological perturbations due to 
Ellis and Bruni. In Sec. 3, we show how the energy conservation equation 
leads naturally to a covector which obeys a conservation equation. In Sec. 4, 
we show explicitly the connection between our approach and the familiar, 
coordinate-based, linear theory. In Sec. 5, we study second order perturba- 
tions. In Sec. 6, we consider other non-linear conserved quantities. In Sec. 7, 
we discuss the gauge invariance of our variables and finally, in Sec. 8, we 
summarize and conclude. 

2 Covariant formalism 

In this section, we briefly review the basic ideas of the covariant approach 
developed by Ellis, Bruni and collaborators (2EJ EE EHj , and based on ear- 
lier works by Hawking jHOj and Ellis |31j . For simplicity, we consider the 
universe filled with a single perfect fluid, although most of what we will say 
applies to any independent (i.e., non- interacting with other fluids) perfect 
fluid, whether it is or not surrounded by other fluids. The case of an im- 
perfect fluid, with non-vanishing heat flow and anisotropic stress, will be 
considered in a separate publication. The fluid we consider can be charac- 
terized by a comoving four-velocity u a = dx a /dr (u a u a = —1), where r is the 
proper time along the flow lines, a proper energy density p and a pressure P. 
The energy-momentum tensor associated to the perfect fluid is given by 

T% = (p + P) u a u b + Pg V (1) 
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To fully characterize the fluid, one needs an equation of state relating P to 
p and, possibly, to other physical quantities if the fluid is not barotropic. 
The spatial projection tensor orthogonal to the fluid velocity u a is defined 

by 

h ab = gab + u a u b , {h\h b c = h a c , h a b u b = 0). (2) 
It is also useful to introduce the familiar decomposition 

V b U a = (T ab + UJab + 7^<S>h ab - U a U b , (3) 

with the (symmetric) shear tensor a ab , and the (antisymmetric) vorticity 
tensor u ab ; the volume expansion 6, is defined by 

e = v a u a , (4) 

while u a is the acceleration, with the dot denoting the covariant derivative 
projected along u a , i.e., ' = u a V ' a . 

The integration of 9 along the fluid world lines with respect to the proper 
time r is 

a= l -JdrQ, (5) 

and can be used to define, for each observer comoving with the fluid, a local 
scale factor S = e a . It follows that 



6 = 3d = 3u a V a a. (6) 

Note that a is defined up to an integration constant for each fluid world 
line. A convenient way to specify this integration constant is to introduce 
some reference hypersurface on which a = 0. The arbitrariness in the choice 
of the integration constant along each world-line is thus translated in the 
arbitrariness in the choice of this reference hypersurface. 

We now wish to define, in a geometrical way, quantities that can be inter- 
preted as perturbations with respect to the FLRW configuration. Since we 
work directly with the clumpy spacetime, an obvious difficulty is how to sep- 
arate homogeneous quantities from their perturbations. A way to get around 
this difficulty was first suggested by Ellis and Bruni |2E]- They introduced 
spatial projections (defined as projections orthogonal to the four-velocity u a ) 
of the covariant derivative of various scalar quantities: of the energy density, 

X a = h a b V bP = D aP , (7) 

of the pressure 

Y a = D a P, (8) 
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and the volume expansion, 

Z a = D a Q. (9) 

Here, following 25^, we also introduce the spatial gradient of the integrated 
expansion, 

W a ee D a a. (10) 

Note that these definitions are purely geometrical and depend only on the 
(physical) four-velocity u a . All these quantities automatically vanish in a 
strictly FLRW spacetime: in this sense we call them perturbations. However, 
they are fully non-perturbative quantities, not restricted to linear order in 
a perturbation expansion. The last quantity, W a , corresponds to the spatial 
gradient of the local number of e- folds a. It was introduced in the non- 
perturbative approach in our previous work 25 , because it plays a crucial 
role in our approach, replacing the familiar spatial curvature perturbation, 
similarly to the separate universe picture |2H E2j- Note that a quantity 
similar to W a has already been used in the linearized covariant theory in 



3 Generalized conserved quantities 

Starting from the above definitions, we will now introduce a quantity that 
naturally satisfies a conservation equation. In [32], it was emphasized that 
the conservation, for adiabatic perturbations and on large scales, of the linear 
curvature perturbation on uniform-density hypersurfaces, usually denoted (, 
introduced in can be derived directly, without resorting to Einstein's 
equations, from the conservation of the energy-momentum tensor, 

v n = o. (ii) 

Here, we will use the same starting point in order to define a non-linear 
generalization of (, and we shall make use of the covariant approach described 
in the previous section. 

Let us consider the projection along u a of the conservation equation ([lip. 
This yields 

p + e(p + p) = o, (12) 

where we remind the reader that the dot represents a covariant derivative 
projected along u a , i.e. ' ee u a V a . If one takes the projected gradient of the 
previous expression one gets 

D a (p) + ( p + P)Z a + Q (X a + Y a ) = 0, (13) 
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where we have used the definitions (J7J), (JBJ) and Q. 

We now concentrate on the first term of Eq. (|T3~j) trying to invert the time 
derivative with the spatial gradient. In order to do so, it will be convenient to 
introduce the Lie derivative along u a of tensors. The corresponding definition 
for a covector, which will be useful below, is (see, e.g., [35J for the general 
definition) 

C-uXa = U C V c Xa + XcV a M C . (14) 

Note that, for scalars, the dot derivation is equivalent to the Lie derivation, 
i.e. p = C u p. Thus, we want to invert the Lie derivation along u a with the 
spatial gradient in the first term of Eq. QlMj) . The two do not commute and 
we have 

D a (p) = h a b V b (u c V c p) = (h a b V b u c )v c p + u c h b V b V cP 

= C u (D a p) - V a u c D cP + (h b V b u c - u b V b h a c ) V c p. (15) 

This expression can be rewritten in a remarkably simple form as 

D a (p) = C u {D a p) - u a p. (16) 

Although we have derived it for the energy density p, this expression of 
course applies to any scalar quantity. In particular, it can be applied to the 
integrated expansion a, thus allowing us to express Z a = D a Q = 3D a a in 
terms of W a = D a ot. This gives 

Z a = 3C u W a - 3ii a a. (17) 

Substituting (|T6|) and (fTTj) in (|T3|) . one gets 

C u X a + 3 (p + P) C u W a + 6 (X a + Y a ) = 0, (18) 

where terms proportional to u a have disappeared because of the continuity 
equation (JT2J). Using the relation 

^ £ -(^)- e < i+ ^ (i9) 

where we have introduced c 2 s = P/p, the non-linear generalization of the 
adiabatic speed of sound, Eq. (fTSj) can be rewritten as an evolution equation 
for the covector 

X a 
C a = W a + a = D a a - ^p, (20) 
3(p + P) p 
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which, as we will see later, can be seen as a generalization of the curvature 
perturbation on uniform density hypersurfaces of the linear theory, usually 
denoted by ( (hence our name ( a ). Equation (fT8|) then yields 

£u( a = - g^fpj ( Y * - • (21) 

On the right hand side, the quantity 

Y a = Y a - c 2 s X a = D a P - -D aP (22) 

P 

can be interpreted as the projected gradient of the non-adiabatic pressure, 
and represents the non-linear generalization of the non-adiabatic pressure 
defined, e.g., in [32] ■ It vanishes for purely adiabatic perturbations, i.e., 
when the pressure P is solely a function of the density p. 

For practical purposes, it is useful to note that, although both ( a and r a 
are defined as linear combinations of spatially projected gradients, one can 
replace them by ordinary gradients, thus having 

a P 
(a = d a a - -d a p, T a = d a P - —d a p. (23) 
P P 

This is indeed possible for any linear combination of the form D a x—{x/v)D a V 
since, for any scalar quantity x, we have 

DaX = d aX + U aX - (24) 

Substituting the definition of T a into Eq. (j21j) we finally obtain 

CuCa = ~3(p~+P) Ta - (25) 

This equation has a form very similar to the conservation equation for ( of 
the linear theory, which will be rederived in the next section. However, in 
our case, the above equation is exact, fully non-perturbative and valid at all 
scales. 

The time derivative that appears in the equation of the linear theory has 
been replaced here by a Lie derivative along u a , which is consistent with the 
fact that we want to describe the conservation properties of a covector along 
the fluid worldline of a comoving observer. This makes Eq. ()25|) remarkably 
simple and improves our original formalism given in j2Hj where we did not 
introduce the Lie derivative. Moreover, the use of the Lie derivative simpli- 
fies the calculations in practice, because one does not need to compute the 
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covariant derivatives. Indeed, the left hand side of Eq. (J23|) can be rewritten, 
starting from the definition of Lie derivative ([25)1 . as 

C-uCa = U C d c Ca + CcO a U C , (26) 

which will be very simple to use in the following. 

A remark is in order here. In our approach, the covector ( a , defined 
as a linear combination of spatial gradients of a and of p, is the crucial 
quantity that enables us to write the conservation equation for perturbations 
in a remarkably simple form. It is interesting to note that the same linear 
combination appears in the work of Rigopoulos and Shellard pE] • Although 
their motivation is, like here, to generalize the conservation equation for 
the linear ( to the non-linear case, their approach differs from ours in two 
respects. First they do not use a covariant formalism but a coordinate based 
approach with an ADM decomposition of the metric. As a consequence, their 
correspondent equivalent of our a is not defined as the integrated expansion 
along fluid worldlines but from the determinant of the spatial metric. In this, 
they follow an older work [2H] which already emphasized the advantage of 
this quantity to generalize the linear £ to non-linear order. Second - and as 
a consequence of the first point - they restrict their analysis to super-Hubble 
scales as in [20J. 

In a similar spirit, i.e., using the ADM decomposition and invoking a 
long-wavelength approximation, the authors of fH] found, not exactly the 
linear combination appearing in JE], but its (spatially) integrated version. 
In contrast with these previous works, as already emphasized, our formulation 
is fully covariant and our non-linear generalization applies to all scales. 



4 Linear theory 

In this section, we compare our approach with the more familiar, coordinate- 
based, linear theory (see e.g., [HE] or [Slj for recent presentations of this 
topic). In order to do so, we first introduce conformal coordinates x M = 
{rj, x 1 } to describe our almost-FLRW spacetime. A prime will denote a partial 
derivative with respect to conformal time rj, ' = d/drj. The background 
spacetime, i.e., at zeroth order, is a FLRW spacetime, endowed with the 
metric 

ds 2 = g^dx^dx" = a 2 —drj 2 + r ) i jdx l dx^ , (27) 

where a = e a is the background scale factor, and filled with a homogeneous 
perfect fluid characterized by the energy density p(rj), the pressure P{rj) and 
the four-velocity 

u» = {1/a, 0,0,0}. (28) 



8 



The perturbed spacetime is described by the perturbed metric 

ds 2 = (g^ + 8g^) dx»dx v , (29) 

with 

5g 00 = -2a 2 A, 5g 0i = a 2 Bi, bg^ = a 2 i^. (30) 
We decompose in the form 

Hij = -2^ 7ij + 2ViVjE + 2V (i Ej ) + 2E%, (31) 

with Efj transverse and traceless, i.e., ViE T ^ = and j^Efj = 0, and Ef 
transverse, i.e., ViE Vl = 0, where Vj denotes the covariant derivative with 
respect to the homogeneous spatial metric 7y (which is also used to lower or 
raise the spatial indices). 

The corresponding matter content is a perfect fluid with perturbed energy 
density and pressure, 

p(r], x*) = p(rj) + 6p(r), x l ), P( V , x l ) = P( V ) + SP( V , x l ) (32) 

and four-velocity 

u li = u li + 5u^, 8u» = {-A/a, v^a}, v k = V k v + v%, (33) 

where vj is transverse, ViV Vt = 0. 

We now wish to write down explicitly the components of our vector ( a 
in this generic coordinate system. We will start from the expression for ( a 
given in ()23|) rather than from its definition (}2T)j) . At zeroth order, i.e., in the 
unperturbed FLRW spacetime, ( a automatically vanishes. At linear order, 
the spatial components are simply 

CP = d^ l \ C (1) -« (1) -|p (1) , (34) 

where we recall that a prime denotes a partial derivative with respect to 
conformal time. 

To compute the component ( , it is useful to note that, for any function 
X, one can write 

D x = u u % diX - u l Uid x, (35) 

where we have used the normalization of u a . Since u l is first order, this 
implies that £g = (although this will not be the case at second order). 

In order to make the link with the usual quantities of the linear theory, 
one needs to reexpress in terms of the metric and matter perturbations. 
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The detailed calculations, up to second order in the perturbations, are given 
in the appendix. Retaining only the first order terms, we obtain 



a ( i) 



-tP + ^V 2 E + ^J dr]V 2 v, (36) 



where the decompositions (|3Tj) and have been used to obtain the second 
line. 

The components of the non-adiabatic term T a = d a P — (P/p)d a p can be 
deduced directly from the components of ( a by substituting P to a. There- 
fore, one finds 

if^ftrw r« = p (1) - j P{1) , (37) 

while the time component vanishes. 

Let us now specialize our equation (J25j) to first order in perturbation 
theory. One can first notice that, at first order, C u (-^ = (j- 1 ' 1 /a. We then 
get 

CP = -^V« (38) 
p+P 1 v ' 

where 7i denotes the conformal Hubble rate, 7i = a' /a = a'. We then note 
that, from Eq. (}36j) . 

a (1) ' = -iP' + ~V 2 (E' + v). (39) 

Consequently, making use of Eqs. and (}3Tj) and getting rid of the <9j 
common to both sides of (1551). one finds 



(_,_ H M) + i vV + „).__^_(p (1) _^ (1) ). (40) 

One recognizes in the parenthesis of the left hand side the familiar quantity 

C = -v> - n^-. (4i) 
p 

Our linearized equation for ( a is thus shown to be completely equivalent to 
the analogous equation obtained directly in the linear theory 

C = -^js5P nad - V (E' + v) , (42) 
p + P 3 
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with 

5Pnad = P(l)-yP(l). (43) 

As one can see, at the linear level, our quantity Q is conserved at all 
scales because our definition of Q automatically includes the Laplacian terms 
that appear on the right hand side of (|4~2|). In other words, whereas our ( a 
coincides with the usual ( on long wavelengths when the spatial gradients 
can be neglected, the two quantities will differ on small scales, our quantity 
being more adapted to the underlying conservation law. 



5 Second order perturbations 

In the following, we decompose any function X in the form 



X(r,, x*) =X(rj) + X (1) ( V , x l ) + -X m (n, x l ) + 



(44) 



where Xn\ and Xr 2 ) represent, respectively, the first and second order per- 
turbations. 

To compute the component Co at second order we use again Eq. (|35|) . Since 
u l is first order, this implies (jp = u l Q^ . Expanding now Q = dia — (a/p)dip 
up to second order, one finds 



Cf ) = di I a {2) - ^-p (2) ] - ^ I a' (1) - ^-p' (1) I d iP{1) 



a 



a' 



P 



P' 



P' 



(45) 



In contrast with the first order expression, the coefficient in front of the gradi- 
ent in the second term on the right hand side cannot be directly absorbed in 
the gradient because it depends on space. One can however do the following 
manipulations: 



(2) 



a' 



Oi I «(2) - — P(2) 



a' 



P 



a 



(1) " fP'iD j 9 iPV) 



a a o / P(l) . _,P(1)P(1) 

dt I a (2) - jP{2) - 2« (1) ^ + 2a -^-^ 



a! 



+2d t ( a[ x) - -p' (1) 



PW 
P 



0; 



"(2) - — P(2) 



(1) p p p 



. / _// - f P \ P(i) 



P 



(46) 
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which can be written in the form 



Cf ) = ^C (2) + y { iAC {1) \ (47) 



with 



. ( o\ a 1 2 . a! . 1 (a! \ 2 . . 

Q ' = "(2) - -P(2) - -"(1) P(i) + 2— p (1) p (1) + - I — I p (1) . (48) 

If one substitutes in the above equation the expression for a up to sec- 
ond order derived in the appendix, keeping only the scalar terms without 
gradients, i.e., which are not negligible on large scales, 

a ~ In a — ip — ip 2 , (49) 

one finds 

T~C 2 T~C 1 ( T~t\ 

C (2) „ _^ (2) _2V; (1) 2 --p (2) + -V; (1) V (1) +2— + P(1) 2 . (50) 

The right hand side can be easily related to the conserved second order 
quantity defined by Malik and Wands in [12], and thus, 



s — SMW smw- 

(See also the discussion in 

Now that we have identified the second order perturbation variable 
we can expand Eq. (|25j) to second order. On making use of Eq. ()26)) to 
reexpress the Lie derivative along u a in terms of coordinate time derivative, 
and retaining only second order terms, we have 

C U (P = - [C {2)/ - 2AC {1) ' + 2 (vWjCV + # V)] • (52) 

CI 

Finally, on making use of Eqs. (JSH), flH|), and (gEJ), and that 9 = 3(l-A)af/a 
up to first order, we can explicitly write the conservation equation (|23j) at 
second order and on all scales: 

(2)/ = ^ r (2)_^_ r (l) C (l)/_ 2 ^ C (l)_ (53) 

p + P p + P J K ' 

The definition of can be read from the expansion of by substituting 
P to a. For adiabatic perturbations, we find that the second order scalar 
variable Q 2 ' is conserved only on large scales, when the last term on the right 
hand side of Eq. (j33J) can be neglected. In principle, one can extend straight- 
forwardly the procedure presented here to higher orders in the perturbation 
expansion. 
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6 Other conserved quantities 



6.1 Comoving curvature perturbation 

We have seen in the previous sections that our ( a can be interpreted as the 
non-linear generalization of the curvature perturbation on uniform energy 
density hypersurfaces of the linear theory. Another useful quantity in the 
linear theory is the so-called curvature perturbation on comoving hypersur- 
faces, usually denoted by 1Z and one can wonder if a non-linear generalization 
can be constructed within our formalism. As we have seen before, what plays 
the role of the spatial curvature in our formalism is the integrated expansion 
a and since our spatial gradients are defined with respect to the comoving 
observers, a candidate that generalizes 7Z is simply 

1Z a = -D a a = -W a . (54) 

We can explicitly check the connection between this quantity and the 
linear comoving curvature perturbation 7Z by noticing that, in the coordinate 
system defined in the appendix, one finds for the spatial components of lZ a , 
at first order in the perturbations, 

Ti? = -<W 1} - am = -dia {1) - a'iyi + B,). (55) 

Therefore, for the scalar part of 7?j , we recover the usual definition of 
the comoving curvature perturbation, 

-rSM = diRM, 1Z W = -H(v + B), (56) 

provided one can replace — cr 1 ) by ij). 

It is relatively easy to derive an evolution equation for 7Z a . As a first 
step, one can rewrite Eq. (JTTj) as 

C u K a = ~ (Qu a + Z a ) . (57) 

If one wishes to write the right hand side of this equation only in terms of 
quantities characterizing the matter, there is some extra work. This consists 
in replacing Z a with the use of the field equations, i.e., Einstein's equations. 
Contracting the identity 

V c V rf M a - Vd^cUa = RabcdU b ', (58) 

where R a bcd is the Riemann curvature tensor, by g ac h ed , one gets, upon using 
the decomposition (J3J and Einstein's field equations [R a b = « {T a b — \Tg a ^j , 
with k = 8nG], the so-called momentum constraint equation (see, e.g., [3*T]). 

2 -Z a - h a b V c {o\ + O + (°a + ^a)ub = -Kh a b T bc u c = 0, (59) 
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where the right-hand side vanishes for the energy- momentum tensor ([T}. 
Moreover, one can write u a as 



4. = (60) 

This is simply the Euler equation, which follows from the projection, or- 
thogonally to u a , of the energy-momentum tensor conservation (jllj) . The 
covector Y a can be reexpressed as 

Y a = T a + c 2 s X a . (61) 

It is also worth noticing the identity that relates our two covectors ( a and 



Finally, substituting the relations (j59H61|) into Eq. (|57|) . we obtain 

b _i_ , , b~ 

(V, -i c 2 X.) - y 



""u'^-a 



3(p + P) 2(p + P) 



{v b + c 2 s X h )- l -h a b V c {a\ + u\). (63) 



This gives a fully non-linear evolution equation for lZ a . As one can see, it is 
much more complicated than the evolution equation for ( a . 

At linear order, considering only scalar perturbations, one finds, from the 
definition of X a , 

X« = p {1) - 3H(p + P)(v + B), (64) 

where Xf^ 1 ' = diX^\ The quantity X^ 1 ^ corresponds to the so-called comov- 
ing density perturbation [2] , known to vanish on large scales because of the 
relativistic Poisson equation. Therefore, since Q 1 ' +1Z^ = X^ /[3(p + P)}, 
—1Z^ coincides with on large scales, and is thus conserved on large 
scales. Note that this also holds at second order, as shown in ITj 



6.2 The conserved variable C a 

In the previous subsection, we introduced the variable lZ a and showed that, 
in the linear limit and on large scales, it reduces to the familiar comoving 
curvature perturbation. This is however not the only non-linear quantity 
that satisfies this property. As shown in [2H], this is also the case for the 
projected gradient of the spatial scalar curvature, 

C a = S 3 D a IC. (65) 
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To be more precise, the quantity K, is denned as 



-±e> 

3 



a' 



to + Kp 



(66) 



with cr 2 = \a a b(y ab and u; 2 = \u) a bUJ ab . When u a b — 0, K corresponds to the 
Ricci scalar curvature of the three-dimensional hypersurfaces orthogonal 
to u a . However, in the general case when oj ah ^ 0, such hypersurfaces cannot 
be defined, according to Froebenius' theorem (see, e.g., [35J). 

As shown in [28J, C a is conserved on large scales in a linearly perturbed flat 
FLRW universe, for adiabatic perturbations. More generally, it is convenient 
to define the quantity [39 j 



C a — Ca 



AaK 
p + P 



(67) 



where K = 0, ±1 characterizes the spatial curvature of the unperturbed 
FLRW spacetime. The variable C a reduces to C a for a spatially flat back- 
ground and turns out to be conserved on large scales for adiabatic linear 
perturbations for all values of K. 

To check the relation between C a and our quantity lZ a at linear order, we 
now work with the perturbed metric of Sec. 0] The first order perturbation 
of K is given by 



4 r 



(V 2 + 3K)tP - HV k (v k + B k ) 



(68) 



It follows from the definition (|65|) that Cq 1 ^ = and, using /C = 6K/ a 



C 



(i) 



■H + Bi)d Q K: + V^/C {1) 
4a {Vi [(V 2 + 310 (V> - W(« + B))] - 3HK(vf + Bj )} . (69) 



On making use of Eq. (|67jl one finds C ( 

CP 



(i) 



and 



4aV 4 1 V 2 [if> - H{v + B)\ + 3K 
Aadi (V 2 H - 3K( 



3(p + P) 



(70) 



where ( and TZ are the standard first order perturbation variables. 



6.3 Conserved number densities 

We can define other types of conserved quantities starting from any scalar 
quantity which obeys a continuity equation. In analogy with the analysis of 
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Lyth and Wands in the context of the linear theory ;4Qj, one can, for exam- 
ple, consider the particle number density in the physical contexts where the 
particle number is conserved. The corresponding non-perturbative continuity 
equation is simply 

V a (nu a ) = 0, (71) 

which yields 

n + Qn = 0. (72) 

Via a derivation similar to that of Sec. 3 for the energy density, it is not 
difficult to show that the spatial projection of (J72*|) can be rewritten as 



£nCi n) = 0, (73) 



where we have defined 



CW = W a + ^ = D a a - ^D a n. (74) 
on n 

We thus see that, for any quantity that satisfies a standard local conservation 
equation, one can construct a covector that embodies a non-linear perturba- 
tion for this quantity, defined as a combination of the spatial gradient of the 
quantity and that of a. This covector then satisfies an exact fully non-linear 
conservation equation. 



7 Non-perturbative approach and gauge in- 
variance 

The approach proposed in this work is completely non-perturbative: it does 
not rely on a perturbative expansion of the variables describing cosmological 
perturbations. This is why the equations derived here encode the non-linear 
evolution of these variables. In this section, we would like to comment about 
the notion of gauge-invariance in the context of this approach. 

Gauge-invariance has become ubiquitous in the studies of linear cosmo- 
logical perturbations. Indeed, using gauge-invariant quantities helps avoiding 
ambiguities by getting rid of unphysical degrees of freedom and makes easier 
the comparison between calculations done in different gauges. However, com- 
pletely fixing a gauge, either physically (i.e., with reference to some specific 
matter) or geometrically, is equivalently acceptable. 

In our case we do not really need to care about gauge-invariance because 
we use tensor quantities that are physically well specified independently of 
any coordinate system. Our approach, however, contains one possible source 
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of ambiguity: our quantity a is defined up to a constant of integration along 
each fluid worldline. Assuming that a is continuous implies that there is some 
arbitrariness in the choice of the a = hypersurface. A similar arbitrariness 
can be found in the separate universe approach 

This being said, it is nevertheless instructive to examine our approach in 
the light of the treatment of gauge-invariance beyond linear perturbations 
which has been discussed thoroughly in [TU]. In the perturbative approach, 
one considers a collection of space-times J\A\ that interpolates continuously, 
as the parameter A goes from to 1, between the background spacetime Aio 
and the "real" universe A4\. A choice of gauge corresponds to a choice of a 
continuous family of diffeomorphisms (p\ : A4q —>■ M.\. The perturbation of 
a tensor field T, which is associated to a continuous family of T A defined on 
each A4\, is defined as 

AT A = iflT - T , (75) 

where tp* x T is the pull-back of T\. The perturbation AT\ is thus a tensor field 
defined on Aio- One can then do a Taylor expansion of the perturbation with 
respect to the parameter A, 

co \k 

AT a = E ^[4 T > ( 76 ) 

where £ is the tangent vector associated with the flow ip\, and one can thus 
define, order by order, the notion of gauge-invariance [TU] . 

In a non-perturbative approach, such as ours, one can define an alter- 
native notion of gauge-invariance. Indeed one needs to consider only the 
background spacetime M.q and the real spacetime M.. A correspondence be- 
tween M.q and M. can be specified via a diffeomorphisms ip : M. Q — > M.. The 
gauge ambiguity results from the arbitrariness in the choice of ip. However, 
now, instead of defining the perturbation of a tensor T on the background 
spacetime .Mo, one can do it on the perturbed spacetime: 

AT = T - <p m T , (77) 

where (p* is the push-forward of T. With this alternative definition for the 
perturbation, it is immediate to see that AT is gauge- invariant, i.e., indepen- 
dent of the choice of <p>, if T vanishes. In this sense, which does not coincide 
with the perturbative gauge-invariance defined in ^U] but which agrees im- 
plicitly with the discussion of [2^1, the quantities defined in our approach are 
gauge- invariant . 
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8 Conclusions 



In this work, we have discussed in detail a new approach to the theory of non- 
linear cosmological perturbations, which was presented briefly in Our 
starting point is the covariant formalism for cosmological perturbations and 
in particular its use of a spatially projected gradient, defined orthogonally 
to the fluid four-velocity Our new formalism is then based on the following 
additional ingredients: 

• We have introduced a, the local number of e- folds, or integrated ex- 
pansion, along each worldline, and its spatially projected gradient. 

• We have introduced ( a , a linear combination of the spatial gradients of 
a and of the energy density p. For linear perturbations and on large 
scales, our definition reduces to the well-known curvature perturbation 
on uniform energy density hypersurfaces (more precisely, its gradient). 

• We have shown that ( a obeys a remarkably simple conservation equa- 
tion, Eq. (|25j) . which is fully non-linear, exact and valid at all scales. 
Here, improving the formulation of our previous work |25j, we describe 
the evolution of ( a by using the Lie derivative with respect to the 
four-velocity u a , instead of the covariant derivative with respect to u a . 
This explicitly shows that ( a is conserved, at all scales and fully non- 
perturbatively, for adiabatic perturbations. Indeed, the only source of 
evolution of ( a is the non-linear generalization of the non-adiabatic pres- 
sure, which describes entropy perturbations and is defined as a linear 
combination of the spatial gradients of p and P. 

We have shown explicitly the connection between our approach and pre- 
vious works, in particular on linear perturbations and on second order per- 
turbations. As an example, our equation (|2*5Jl is so simple that it is straight- 
forward to write down the conservation equation for the components of ( a at 
second order and at all scales. One of the conclusions of this work is that the 
local number of e-folds and the related combination ( a appear to be the most 
natural quantities to express the conservation law for adiabatic cosmological 
perturbations, in its fully non-linear and all-scale form. 

As a final remark, let us stress that we did not use Einstein's equations 
to derive the conservation equation for ( a but only the conservation of the 
energy-momentum tensor. Therefore, our equation (f23j) applies to any grav- 
ity theory, as long as the energy-momentum tensor is conserved. 
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A Appendix: Perturbed quantities up to sec- 
ond order 

We consider a perturbed metric with 

#00 = -a 2 (l + 2 A), g 0i = a 2 Bi, g tj = a 2 (7^ + H^) , (78) 

where A, B t and are the perturbations. We decompose in the form 

Hij = -2^ + 2V i V j £ + 2V (l Ej } + 2Ej, (79) 

with Efj transverse and traceless (TT), i.e. ViE Tl ^ = and Y^Efj = 0, and 
Ej is transverse, i.e. ViE v% = 0. The components of the inverse metric g ab 
are given, up to second order, by the following expressions: 



9°° 

a 



= ^ (-1 + 2A- 4A 2 + BiB*) , g 0i = 4 (B i - 2AB 1 - H ik B k ) , 
a 2 v J a 2 v ' 

= — (yi - H ij - B i B j + H\H kj ) , (80) 
a 2 v J 

where the spatial indices are raised or lowered via the background comoving 
spatial metric 7^. One has to imagine that each perturbation quantity can 
be expanded up to second order, e.g., A = A^ + \A^ 2 \ 

The time component u° of the four-velocity can be derived from the nor- 
malization condition g a bU a u b — — 1, and one finds at second order 

u° = \ (l - A + \A 2 + l -v k v k + B k v k ) , (81) 

where we have introduced v k = au k . 

We now have all the ingredients to compute the expansion = V a u a up 
to second order. One finds 

6 = l{(l-A)(^n + lH\^-^H ij H^ + l[(v k + B k )(v k + B k )]' 
+m (I a 2 + \ VkV k + BkV kj + VkV k + v k Vk ^ + 1 ^ j ) (g2) 
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where H\ = = —Gip + 2V 2 E. This formula must be compared to the 

expression of 9 as a function of the integrated expansion a, = 3d = u a V a a, 
which reads up to second order: 



e = ? 

a 



3 1 

1 - A + -A 2 + -v k v k + B k v h ) a' + v k V k a 



(83) 



The order by order comparison finally yields 

i ri„„, i i 



a' 



n + - {-H\' - -HvH* + - {(v k + B k ) (v k + B k ] 



+ (l A)V k v k + v k V k (A+^H%)\ /■'' V/.n'". (SI) 
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